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Abstract 

The proton spin puzzle inspired by the EMC experiment and its present status 
are closely examined. Recent experimental progress is reviewed. Various factorization 
schemes due to the ambiguity arising from the axial anomaly are discussed. Some 
misconceptions in the literature about the MS factorization scheme are clarified. It is 
stressed that the polarized nucleon structure function g\ (x) is independent of the factor- 
ization scheme chosen in defining the quark spin density. Consequently, the anomalous 
gluon and sea-quark interpretations for the deviation of the observed first moment of 
g±(x) from the Ellis-Jaffe sum rule are equivalent. While it is well known that the total 
quark spin in the chiral-invariant (CI) factorization scheme (e.g. the improved parton 
model) can be made to be close to the quark model expectation provided that the gluon 
spin is positive and large enough, it is much less known that, contrary to the gauge- 
invariant scheme (e.g. the MS scheme), the quark orbital angular momentum in the 
CI scheme deviates even farther from the relativistic quark model prediction. Recent 
developments in the NLO analysis of polarized DIS data, orbital angular momentum 
and lattice calculations of the proton spin content are briefly sketched. 
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1 Introduction 



Experimentally, the polarized structure functions g\ and #2 are determined by measuring 
two asymmetries: 

da^ — da^ da^ — do^~* , . 

11 = dan + da^ ' X = da^ + da^ ' { ' } 

where da^ (da^) is the differential cross section for the longitudinal lepton spin parallel 
(antiparallel) to the longitudinal nucleon spin, and da^ (dcr>~~*) is the differential cross 
section for the lepton spin antiparallel (parallel) to the lepton momentum and nucleon spin 
direction transverse to the lepton momentum and towards the direction of the scattered 
lepton. From the parton-model or from the OPE approach, the first moment of the polarized 
proton structure function 



IW)= r 9 i{x,Q 2 )dx, (1.2) 
Jo 



can be related to the combinations of the quark spin components via 

r?= Iy,4 a q(Q 2 ) = lE4(p,s\qwq\ P ,s)s», (1.3) 

where Aq represents the net helicity of the quark flavor q along the direction of the proton 
spin in the infinite momentum frame: 

f 1 f 1 r 1 

Aq = / Aq(x)dx = / q\x) + <f{x) — q^{x) — q^{x) dx. (1.4) 



At energies (Q 2 ) ~ 10 GeV 2 or smaller, only three light flavors are relevant: 

IW) = ~ (^A M (Q 2 ) + l -Ad{Q 2 ) + l -As{Q 2 )) . (1.5) 

Other information on the quark polarization is available from the low-energy nucleon axial 
coupling constants g\ and g\: 

g\{Q 2 ) = (p,s\u^ 5 u-d^ 5 d\p,s)s^ Au(Q 2 )-Ad(Q 2 ), (1.6) 
g s A (Q 2 ) ee (p, s\u lfl j 5 u + djtf B d - 2s~?^ 5 s\p, s)s» = Au(Q 2 ) + Ad(Q 2 ) - 2As(Q 2 ). 

Since there is no anomalous dimension associated with the axial- vector currents A?, and Af., 
the non-singlet couplings g\ and g\ do not evolve with Q 2 and hence can be determined at 
q 2 = from low-energy neutron and hyperon beta decays. Under SU(3)-flavor symmetry, 
the non-singlet couplings are related to the SU(3) parameters F and D by 

g 3 A =F + D, g 8 A =3F-D. (1.7) 

We use the updated coupling g\ = 1.2670 ± 0.0035 [JTJ] and the values @ 

F = 0.463 ± 0.008 , D = 0.804 ± 0.008 , F/D = 0.576 ± 0.016 (1.8) 
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to obtain g\ = 0.585 ± 0.025 . 

Prior to the EMC measurement of polarized structure functions, a prediction for 
was made based on the assumption that the strange sea in the nucleon is unpolarized, i.e., 



As = 0. It follows from ( |1.5|) and ( |1.6| ) that 



r?(Q 2 )= ^9 3 a + ^9% (1.9) 

This is the Ellis- Jaffe sum rule 0: rf = 0.185 ± 0.003 in the absence of QCD corrections 
and equals to 0.171 ± 0.006 at Q 2 = 10 GeV 2 to leading-order corrections. The 1987 EMC 
experiment Q then came to a surprise. The result published in 1988 and later indicated that 
rf = 0.126 ± 0.018, substantially lower than the expectation from the Ellis- Jaffe conjecture. 
From the EMC measurement of rf, we obtain 



Am = 0.77 ±0.06, Ad = -0.49 ± 0.06 , As = -0.15 ± 0.06 , (1.10) 



and 



9a(Q 2 ) = (p, s|«7/i75U + d^ 5 d + s-f^ 5 s\p,s)s^ 

= Au{Q 2 ) + Ad(Q 2 ) + As(Q 2 ) = AS(Q 2 ) = 0.14 ± 0.18 (1.11) 

at Q 2 = 10.7 GeV 2 . The EMC results exhibit two surprising features: The strange-sea 
polarization is sizeable and negative, and the total contribution of quark helicities to the 
proton spin is small and consistent with zero. This is sometimes referred to as the "proton 
spin crisis". 

The so-called "proton spin crisis" is not pertinent since the proton helicity content ex- 
plored in the DIS experiment is, strictly speaking, defined in the infinite momentum frame 
in terms of QCD current quarks and gluons, whereas the spin structure of the proton in 
the proton rest frame is referred to the constituent quarks. That is, the quark helicity Aq 
defined in the infinite momentum frame is generally not the same as the constituent quark 
spin component in the proton rest frame, just like that it is not sensible to compare apple 
with orange. What trigged by the EMC experiment is the "proton helicity decomposition 
puzzle" rather than the "proton spin crisis" (for a review, see || ^, |7], ||, |, [10], |TT)). The 
non-relativistic SU(6) constituent quark model predicts AS' = AU + AD = 1 (AU denoting 
the constituent up quark spin and likewise for AD), but its prediction for the axial- vector 
coupling constant g\ = | is too large compared to the measured value of 1.2670 ± 0.0035 [|I|]. 
In the relativistic quark model, the proton is no longer a low-lying S*-wave state since the 
quark orbital angular momentum is nonvanishing due to the presence of quark transverse 
momentum in the lower component of the Dirac spinor. Realistic models e.g. the cloudy 
bag model |12[], predict AS' w 0.60; that is, about 40% of the proton spin is carried by 
the orbital angular momentum of the constituent quarks. In the parton picture, the naive 
expectation of AS, which is equal to g\ = 0.59 under the assumption of vanishing sea po- 
larization, is very close to the relativistic quark model's prediction of AS'. One of the main 
theoretical problems is that hard gluons cannot induce sea polarization perturbatively for 
massless quarks due to helicity conservation. Hence, it is difficult to accommodate a large 
strange-sea polarization in the naive parton model. 
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2 Experimental Progress 



Before 1993 it took 5 years on the average to carry out a new polarized DIS experiment 
(see Table I). This situation was dramatically changed after 1993. Many new experiments 
measuring the nucleon and deuteron spin-dependent structure functions became available. 
The experimental progress is certainly quite remarkable in the past years. 

Since experimental measurements only cover a limited kinematic range, an extrapolation 
to unmeasured x — > and x — > 1 regions is necessary. At small x, a Regge behavior g±(x) oc 
x a{o) j g conventionally assumed in earlier experimental analyses. The improvement by the 
new measurements is two-folds: First, the small x region has been pushed down to the order 
10~ 3 in SMC experiments (see Table I). Second, an extrapolation to the unmeasured small 
x region is done by performing a NLO QCD fit rather than by assuming Regge behaviour. 
The uncertainties of 17 which mostly arise from the small x extrapolation are substantially 
reduced. From Table I it is also clear that the EMC experiment has been confirmed by all 
successive polarized DIS measurements. 

Table I. Experiments on the polarized structure functions g\(x,Q 2 ), gi(x,Q 2 ) and gf(x,Q 2 ). 



ptargci , ! j j ; 



Experiment Target 
E80/E130 @,|14| p 
E142 [pi 
E143 [16| 



Q (GeV ) range x range 



\Q 2 ) 



E154 |7 
E155 |i 



n 

p,d 
n 

p,d 



KQ 2 
1 < 
1 < 



< 10 
2 < 10 
2 < 10 



1< Q 2 < 17 
1< Q 2 < 17 



O.K x < 0.7 
0.03 < x < 0.6 
0.03 < x < 0.8 

0.014 < x < 0.7 
0.01 < x < 0.9 



17(10) = 0.17 ±0.05* 
17(2) = -0.031 ± 0.006 ± 0.009 
rf (3) = 0.132 ± 0.003 ± 0.009 
Tf (3) = 0.047 ± 0.003 ± 0.006 
17(5) = -0.041 ± 0.004 ± 0.006 
Tf (5) = 0.0266 ± 0.0025 ± 0.0071 



EMC U 
SMC PI 



P 

p,d 



1< Q 2 < 200 
1< Q 2 < 60 



0.01 < x < 0.7 
0.003 < x < 0.7 



(10.7) = 0.126 ± 0.010 ± 0.015** 
rf (10) = 0.120 ± 0.005 ± 0.006 ± 0.014 
17(10) 

r? (io) 



-0.078 ± 0.013 ± 0.008 ± 0.014 
0.019 ±0.006 ±0.003 ±0.013 



HERMES [20] 



p,n,d 



1<Q 2 < 10 0.023 < x < 0.6 17(3) = -0.037 ± 0.013 ± O.OOSt 



* Obtained by assuming a Regge behavior A\ oc x 114 for small x. 

** Combined result of E80, E130 and EMC data. The EMC data alone give r* = 0.123 ± 0.013 ± 0.019 . 
t gf(x)dx = 0.122 ± 0.003 ± 0.010 is obtained by HERMES for the proton target. 



Comparing to the original EMC measurement, the statistic and systematic errors of 
the combined world average for are substantially reduced. The result is = (0.12 ~ 
0.13) ± 0.007 at Q 2 = (5 ~ 10) GeV 2 . Consequently, AS = (0.20 ~ 0.30) ± 0.04. For 
example, AS = 0.25 ± 0.04 leads to 



Am = 0.81 ±0.01, Ad 



-0.45 ±0.01, As 



-0.11 ±0.01. 



We will employ ( |2.1| ) as the benchmarked values for Ag in ensuing discussions. 
The Bjorken sum rule evaluated up to for three light flavors is [^TJ 



r?(Q 2 ) - r?(Q 2 ) = i 9 -± 

6 g v 



1 - 



43 fa s (Q 2 )\ 2 



12 



20.22 



71 



71 



(2.1) 



(2.2) 
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A serious test of the Bjorken sum rule, which is a rigorous consequence of QCD, became 
possible since 1993. The current experimental results are 



E143 [16] : rf-r? = 0.164 ±0.021, 

SMC || : r? - r? = 0.174+°;°^ ; (2.3) 
at Q 2 = 5 GeV 2 , to be compared with the prediction 

rf-r™ = o.i8i ±0.003 (2.4) 

at the same energies. Therefore, the Bjorken sum rule has been confirmed by data to an 
accuracy of 10% level. 

The quark polarization Aq is usually determined from the inclusive data of g\ by assuming 
SU(3) flavor symmetry. Moreover, inclusive DIS determines the sum of polarized quark 
and antiquark distributions, but not the valence and sea quark spin distributions. Semi- 
inclusive polarized experiments in principle allow the determination of Aq for each flavor 
and disentangle valence and sea polarizations separately p2[ . Hence, SU(3) flavor symmetry 



can be tested by comparing the measured first moments of the flavor distributions to the 
SU(3) predictions P3"| . Semi-inclusive data are available from SMC |2l| and HERMES 



(see Table II). In order to present the experimental results, we digress for a moment 
to adopt a different definition for quark spin densities here: Aq s (x) = ql(x) — q^(x) and 
Aq(x) = q\x) — q^(x), where Aq s (x) = Aq(x) — Aq v (x) is the sea spin distribution for the 
quark flavor q. The SMC analysis is based on the assumption of SU(3) flavor symmetric sea: 
A«(a;) = Ad(x) = As(x) = Au s (x) = Ad s (x) = As(x), while the HERMES results shown 
in Table II rely on the assumption of flavor independent polarization: 

Au s (x) _ Ad s (x) _ As(x) _ Au(x) _ Ad(x) _ As(x) 
u s {x) d s (x) s(x) u{x) d(x) s(x) 

Note that the ansatzQ of Aq(x) = Aq s (x) has to be made in both experiments in order 
to extract the valence quark polarization Aq v from the measurement of Aq and Aq, i.e., 
Aq v (x) = Aq(x) — Aq(x). However, one caveat has to be mentioned: A priori the antiqaurk 
spin Aq can be different from the sea polarization Aq s if they are not produced from gluons. 



For example, antiquarks are not polarized in the model of Under the assumption of 
SU(3) symmetric sea polarization, we are led to the predictions g\ = Au v — Ad v , g\ = 
Au v + Ad v and hence Au v = 2F = 0.93 ± 0.02 and Ad v = F-D = -0.34 ± 0.02. Note that 
the valence polarization Aq v should be scale independent. 

The measurement of the gluon spin AG by all possible means is very important both 
theoretically and experimentally (see M for various processes sensitive to the gluon spin 
distributions). A global fit to the present inclusive DIS data of g\{x) cannot even fix the sign 
of AG decisively (see Sec. 3.3), not mentioning its magnitude. One way of measuring AG(x) 
directly is via the photon gluon fusion process occurred in the semi-inclusive DIS reaction. A 

lr rhis ansatz is generally not fulfilled by the assumption of flavor independent polarization made by 
HERMES & 



5 



Table II. The SMC 24] and HERMES [23] results for the first moments of valence and sea spin 
distributions. 





SMC (Q 2 = lOGeV 2 ) 


HERMES (Q 2 = 2.5 GeV 2 ) 


Au v 


0.77 ±0.10 ±0.08 


0.57 ±0.05 ±0.08 


Ad v 


-0.52 ±0.14 ±0.09 


-0.22 ±0.11 ±0.13 


Au 


0.01 ±0.04 ±0.03 


-0.01 ±0.02 ±0.03 


Ad 


0.01 ±0.04 ±0.03 


-0.01 ±0.03 ±0.04 


xAu v 


0.155 ±0.017 ±0.010 


0.13 ±0.01 ±0.01 


xAd v 


-0.056 ±0.026 ±0.011 


-0.02 ±0.02 ±0.02 



recent HERMES |26] measurement of the longitudinal spin asymmetry in photoproduction 
of pairs of hadrons with high transverse momentum indicates that (AG(x)/G(x)) = 0.41 ± 
0.18±0.03 at (x) = 0.17 to LO QCD. Hence AG(x) is found to be positive in the intermediate 
x region. 

3 Theoretical Progress 

In my opinion there are four important progresses in theory: 

• The role played by the gluon to the first moment of gi is clarified. There are two 
extreme factorization schemes of interest. 

• First-principles calculations of the quark spin and orbital angular momentum by lattice 
QCD became available. 

• A complete and consistent NLO analysis of g\ data became possible. 

• Evoluation and gauge dependence of the quark orbital angular momentum are explored. 

3.1 Anomalous gluon and sea quark interpretations 
3.1.1 Anomalous gluon interpretation 

We see from Sec. II that the polarized DIS data indicate that the fraction of the proton 
spin carried by the light quarks inside the proton is AS = (0.20 ~ 0.30) and the strange- 
quark polarization is As « —0.10 at Q 2 = (5 ~ 10) GeV 2 . The question is what kind 
of mechanism can generate a sizeable and negative sea polarization. It is difficult, if not 
impossible, to accommodate a large As in the naive parton model because massless quarks 
and antiquarks created from gluons have opposite helicities owing to helicity conservation. 
This implies that sea polarization for massless quarks cannot be induced perturbatively from 
hard gluons, irrespective of gluon polarization. It is also unlikely that the observed As comes 
solely from nonperturbative effects or from chiral-symmetry breaking due to nonvanishing 
quark masses. 
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As an attempt to understand the polarized DIS data, we consider QCD corrections to 
the polarized proton structure function g[(x). To the next-to-leading order (NLO) of a s , the 
expression for g±(x) is 

9i(x,Q 2 ) = ~Y,4[ AG «( x > a s)® A Q( x >Q 2 )+ AG G( x , a s)®&G(x,Q 2 ) 

= ^E e ?[ A ^ (0) (^Q 2 )+A? (1) (^Q 2 )+A^(^Q 2 ) (3-1) 
+ AC« (x, a.) ® Aq<® (x, Q 2 ) + AC^ (x, a.) ® AG(x, Q 2 ) + ■■•], 

where uses of AC^(x) = 5(1 - x) and ACf {x) = have been made, ® denotes the 
convolution 




and AC q , ACg are short-distance quark and gluon coefficient functions, respectively. More 
specifically, ACq arises from the hard part of the polarized photon-gluon cross section, 
while AC^ from the short-distance part of the photon-quark cross section. Contrary to the 
coefficient functions, Aqf(x) and Aq^(x) come from the soft part of polarized photon-gluon 
and photon-quark scatterings, respectively. Explicitly, they are given by 

Aq {1) (x,Q 2 ) = A^ q (x)®Aq^(x,Q 2 ), Aqf(x,Q 2 ) = A^] G {x) ® AG(x, Q 2 ), (3.3) 

where A<f>j/i(x) is the polarized distribution of parton j in parton i. Diagrammatically, 
and A<jy7f G are depicted in Fig. 1. 



.(') 





lis 




Figure 1: Diagrams for the quark spin distributions inside the parton: A<f)^) g and A<p^)q. 

The photon-gluon scattering box diagram is ultraviolet finite but it depends on the choice 
of infrared and collinear regulators. However, the hard part of the box diagram should be 
soft-regulator independent. In the improved parton model, this is done by introducing a 
factorization scale /if act so that the region k\ > /i 2 act contributes to the hard photon-gluon 
cross section. The results are (see e.g. M) 



ACn*,Q 2 ,/4 ct ) C i 

f 1 dxAC ( a\x) C i 
Jo 



(2x 



In 



A'fact 



In 



x 



X 



(3.4) 
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where the reason for introducing the subscript "CI" will become clear below. Hence, 



dxgi(x, Q z 



7T 



E 



Aq t {Q* 



J CI 



2tt 



-AG{Q Z 



(3.5) 



The (1 — 24 ) term in Eq. (|3.5| ) comes from the QCD loop correction, while the a s AG term 
arises from the box diagram of photon-gluon scattering. If the gluon polarization inside 
the proton is positive, a partial cancellation between Agci and f^AG will explain why the 
observed is smaller than what naively expected from the Ellis- Jaffe sum rule. Note that 
unlike the usual QCD corrections, the QCD effect due to photon-gluon scattering is very 
special: The term a s AG is conserved to the leading-order QCD evolution; that is, AG grows 
with InQ 2 , whereas a s is inversely proportional to InQ 2 . 
As a consequence, Eq. (J2.1p is modified to 



a-.. 



Au cl - --AG 

Z7T 



Add 



—AG 

2tt 

—AG 

2tt 



0.81 ±0.01, 
-0.45 ±0.01 
-0.11 ±0.01 



(3.6) 



and 



9°A 



(Am + Ad + As)a 



3a 6 
~2t^ 



AG = 0.25 ±0.04 



(3.7) 



at Q 2 = 5 ~ lOGeV 2 . Eqs. and Q imply that in the presence of anomalous 

gluon contributions, AEci is not necessarily small while Asqi is not necessarily large. In 
the absence of sea polarization and in the framework of perturbative QCD, it is easily seen 
that AG ~ As(2n/a s ) ~ 2.5 at Q 2 = 10 GeV 2 and AS CI ~ 0.58. It thus provides a nice 
and simple solution to the proton spin puzzle: This improved parton picture is reconciled, 
to a large degree, with the constituent quark model and yet explains the suppression of Tf, 
provided that AG is positive and large enough. This anomalous gluon interpretation of the 
observed T^, as first proposed in J27|, |28], (see also |3(J), looks appealing and has become 
a popular explanation since 1988. 

Some historical remarks are in order: 

• Long before the EMC experiment, there already existed three theoretical calculations 
of the photon-gluon box diagram. Kodaira [JT| was the first one to compute the 
moments of the structure functions g\p. Since he worked in the OPE framework, 
there is no decomposition of g° A in terms of quark and gluon spin components. The 



anomalous gluonic contribution to T\ was first put forward by Lam and Li |32| in 1982. 
A calculation of the gluonic coefficient function using the dimensional regularization 



was first made by Ratcliffe [ 33 1 



The original results for the photon-gluon scattering cross section obtained by [p7 , |28| , p9 | 
are not perturbative QCD reliable as they depend on the choice of soft regulators. The 
first moment of ACq\x) is equal to — a s /(27r) in |29[] but vanishes in the work of 
28] , [3~3]| . After the soft part below the factorization scale /if act is removed, the gluon 
coefficient function is given by Eq. (|3.4j) which is soft-cutoff independent. 
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3.1.2 Sea quark interpretation 



According to the operator product expansion (OPE), the moments of structure functions 
can be expressed in terms of hard coefficients which are calculable by perturbative QCD 
and forward matrix elements of local gauge-invariant operators which are nonperturbative 
in nature: 

l 1 dxx n - x g{x) =Y,CUQ 2 )(N\OU0)\N)- (3-8) 

Jo T 

It turns out that there is no gauge-invariant twist-2, spin-1 local gluonic operator for the 
first moment of gi(x) (see e.g. |3"3||). Here we face a dilemma here: On the one hand, the 
anomalous gluon interpretation sounds more attractive and is able to reconcile to a large 
degree with the conventional quark model; on the other hand, the sea-quark interpretation 
of relies on a more solid theory of the OPE. In fact, these two popular explanations for 
the g\ data have been under hot debate over many years before 1995. Though the OPE 
approach is model independent, it faces the questions of what is the deep reason for the 
absence of gluonic contributions to Y\ and how are we going to understand a large and 
negative strange-quark polarization ? 



3.1.3 Factorization scheme dependence 

In spite of much controversy on the aforementioned issue, this dispute was actually resolved 



almost a decade ago |35|, |36| . The key point is that a different interpretation for Y\ corre- 
sponds to a different factorization definition for the quark spin density and the hard photon- 
gluon cross section. The choice of the "ultraviolet" cutoff for soft contributions specifies the 
factorization convention.^ More specifically, since Acfi^ in Eq. (3.1) is ultraviolet divergent, 
it is clear that, just like the case of unpolarized deep inelastic scattering, the coefficient 
functions AC q and ACq depend on how the parton spin distributions Acjy^l are defined, or 

how the ultraviolet regulator is specified on A<p^ . That is, the ambiguities in defining Atfity 

and A<PqJ G are reflected on the ambiguities in extracting AC^ and ACq . Consequently, 
the decomposition of the photon-gluon and photon-quark cross sections into the hard and 
soft parts depends on the choice of the factorization scheme and the factorization scale /z. 
Of course, the physical quantity g\{x) is independent of the factorization prescription (for a 
review on the issue of factorization, see ||). 

However, the situation for the polarized DIS case is more complicated: In addition to all 
the ambiguities that spin-averaged parton distributions have, the parton spin densities are 
subject to two extra ambiguities, namely, the axial anomaly and the definition of 75 in n 
dimension. It is well known that the polarized triangle diagram for A<p q l J G (see Fig. 1) has an 
axial anomaly. There are two extreme ultraviolet regulators of interest. One of them, which 

2 It is misleading to identify the regularization scheme for soft divergences, e.g. the off-shell scheme, with 
the factorization scheme; the former is merely employed to get rid of infrared and collinear divergences 
appearing in the calculation of partonic cross sections. However, the hard gluon coefficient functions are 
soft-cutoff independent, and they depend on the ultraviolet regulator on the triangle diagram for defining 
the polarized quark distribution inside the gluon. In other words, it is the choice of ultraviolet regulator 
rather than the soft cutoff that specifies the factorization scheme. 



9 



we refer to as the chiral-invariant (CI) factorization scheme (sometimes called as the "jet 
scheme" [0, the "parton model scheme" [11| or the u kj_ cut-off scheme"), respects chiral 



symmetry and gauge invariance but not the axial anomaly This corresponds to a direct 
brute-force cutoff ~/ion the k± integration in the triangle diagram ( i.e. k\ < fi 2 ) with k± 
being the quark transverse momentum perpendicular to the virtual photon direction. Since 
the gluonic anomaly is manifested at k\ — > oo, it is evident that the spin-dependent quark 
distribution [i.e. Aq^(x)] in the CI factorization scheme is anomaly-free. Note that this is 
the fc^-factorization scheme employed in the usual improved parton model |27|, |28], p9f . 



The other ultraviolet cutoff on the triangle diagram of Fig. 1, as employed in the approach 
of the OPE, is chosen to satisfy gauge symmetry and the gluonic anomaly. As a result, chiral 
symmetry is broken in this gauge-invariant (GI) factorization scheme and the sea polarization 
is perturbatively induced from hard gluons via the anomaly. This perturbative mechanism 
for sea quark polarization is independent of the light quark masses. A straightforward 
calculation gives pBL p9| 



Atf; G (x) GI = A<f>% (x)a ~ -(1 - x), (3.9) 



a 



71 



where the term 24 (1 — x) originates from the QCD anomaly arising from the region where 
k\ — > oo. Two remarks are in order. First, this term has been erroneously claimed]] in some 



literature p], [IT], 0, R7| to be a soft term coming from kj_ ~ m . Second, although the 



quark spin distribution inside the gluon A<f>g) G (x) cannot be reliably calculated by pertur- 
bative QCD, its difference in GI and CI schemes is trustworthy in QCD. Since the polarized 
valence quark distributions are fc^-factorization scheme independent, the total quark spin 
distributions in GI and CI schemes are related via Eqs. (|3.3|) and (|3.9|) to be 



Aq(x, Q 2 ) Gl = Aq(x, Q 2 )ci - (1 - x) ® AG(x, Q 2 ). (3.10) 

71 

For a derivation of this important result based on a different approach, namely, the nonlocal 
light-ray operator technique, see Muller and Teryaev @ . 

The axial anomaly in the box diagram for polarized photon-gluon scattering also occurs 
at k\ — > oo, more precisely, at k\ = [(1 — x)/Ax]Q 2 with x — > 0. It is natural to expect that 
the axial anomaly resides in the gluon coefficient function ACq^ in the CI scheme, whereas 



3 Some misconceptions in the literature about the MS scheme have to be clarified. It has been argued pH 
that the GI scheme is pathologic and inappropriate since ACq (x)gi, which is "hard" by definition, contains 



an unwanted "soft" term proportional to (1 — x) [see Eq. ( 3.12 )]. The cross section of the photon-gluon 
box diagram contains a term proportional to (1 — x) if the infrared and collinear divergences are regulated 
by the quark mass or by the dimensional regulator. If the ultraviolet regulator for the triangle diagram is 
chirality-preserving, the (1 — x) term, which arises from the region where k\ ~ m^, does not contribute to the 
hard gluon coefficient, as it should be. However, if the ultraviolet regulator preserves the axial anomaly and 
gauge invariance, for example, the MS regulator, chirality will be broken and the axial anomaly is absorbed 
in the quark spin density. It turns out that the effect of the axial anomaly, which is manifested at k\ ~ Mf act 
(the factorization scale), has the x dependence of the (1 — x) form. This explains why ACq in the MS scheme 
has a term proportional to (1 — x) because the axial-anomaly effect must be subtracted from the previous 
gluon coefficient function computed in the chiral-invariant scheme. As pointed out in [j39[ || and again in 
|12j, the (1 — x) term in the gluonic coefficient function in the MS scheme is purely "hard" , contrary to what 
has been claimed previously. 
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its effect in the GI scheme is shifted to the quark spin density. Since ACq\x) is the hard 
part of the polarized photon-gluon cross section, which is sometimes denoted by gi(x), the 
polarized structure function of the gluon target, we have 

ACg>(*) = g?(x)-A<f>%(x). (3.11) 

It follows from Eqs. and ( gig ) that 

A4 1} (x) GI = a4 1} (x) C i + -(1 - x). (3.12) 

7T 

The first moments of ACq{x), J2 q Aq(x) and g\{x) are given by 



C dxAC%\x) Gl = 0, f 1 dxAC^\x) cl = 
Jo Jo Zn 

A£ GI (Q 2 ) = A£ CI (Q 2 ) - nfa ; m AG(Q 2 ), (3.13) 

Z7T 



and 



It = ^E^(Agci(g 2 )-^AG(g 2 )) 

= ^E^Ag GI (g 2 ), (3.14) 

where we have neglected contributions to g\ from A(p^) q and AC^\ Note that ASqi(Q 2 ) 
is equivalent to the singlet axial charge (p,s\J^\p,s). The well-known results ( |3.10| - |3~T4"|) 
indicate that receives anomalous gluon contributions in the CI factorization scheme (e.g. 
the improved parton model), whereas hard gluons do not play any role in in the GI 
scheme such as the OPE approach. From ( |3.14| ) it is evident that the sea quark or anomalous 
gluon interpretation for the suppression of observed experimentally is simply a matter of 
convention 

The MS scheme is the most common one chosen in the GI factorization convention. The 
so-called Adler-Bardeen (AB) factorization scheme often adopted in the literature is obtained 
from the GI scheme by adding the ^-independent term — a s /(2ir) to AC^ 1 via 

cy 

AC G (x) AB = AC G (x) Gl - -i (3.15) 

while 

Aq(x, Q 2 ) AB = Aq(x, Q 2 ) GI + t — AG(y, Q 2 ). (3.16) 

Z7T Jx y 

In general, one can define a family of schemes labelled by the parameter a pB 



Aq(x) a = Aq GI (x) + [(2x - l)(a - 1) + 2(1 - x)} <g> AG(x), (3.17) 

2lY 



which satisfy the relation 



AS a = A£ GI + ^AG, (3.18) 

Z7T 
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to the first moment, but differ in their higher moments. The AB scheme corresponds to a = 2 
and the CI scheme to a = 1. Since the x dependence of the axial-anomaly contribution in the 
quark sector is fixed to be (1— x), it is obvious that all the schemes in this family including the 
AB scheme cannot consistently put all hard anomaly effects into gluonic coefficient functions 
unless a = 1 , contrary to the original claim made in f41|. 



Finally, it should be stressed that the quark coefficient function AC^(x) in the dimen- 
sional regularization scheme is subject to another ambiguity, namely, the definition of 75 in n 
dimension used to specify the ultraviolet cutoff on ^4> q / q (see Fig. 1). For example, AC q 1 \x) 
calculated in the 75 prescription of 't Hooft and Veltman, Breitenlohner and Maison (HVBM) 
is different from that computed in the dimension reduction scheme. The result 

ACf (z) = Cf (*)-^(l + x) (3.19) 

usually seen in the literature is obtained in the HVBM scheme, where C q (x) is the unpolarized 
quark coefficient function. Of course, the quantity Aq^(x) + AC g (x) <8> Aq^(x) and hence 
g\(x) is independent of the definition of 75 in dimensional regularization. 

3.1.4 A brief summary 

In addition to all the ambiguities that spin-averaged parton distributions have, the parton 
spin densities are subject to two extra ambiguities, namely, the axial anomaly and the 
definition of 75 in n dimension. There are two extreme cases of interest: Either the hard 
axial anomaly is manifested in the matrix elements of the quark current (GI scheme) or it is 
absorbed in the gluonic coefficient function so that the quark matrix element is anomaly-free 
(CI scheme). It should be stressed that in the so-called AB scheme, not all hard anomaly 
effects are lumped into ACc and hence the corresponding quark matrix element is not 
anomaly-free. Of course, it appears that the CI scheme is close to the intuitive parton 
picture as the quark spin distribution which does not contain hard contributions from the 
anomaly is scale independent. (The price to be paid is that Ag CI cannot be expressed as 
the matrix element of a local gauge-invariant operator and hence it is difficult to compute 
by lattice QCD. Also, a priori there is no reason that AS should have a simple quark 
interpretation |p5| .) Nevertheless, physically and mathematically GI and CI schemes are 



equivalent. 

Two remarks are in order, (i) It is worth emphasizing that although the suppression 
of T{ can be accommodated by anomalous gluon and/or sea quark contributions, no quan- 
titative prediction of can be made. An attempt of explaining the smallness of g\ has 
been made in the large- iV c approach (ii) So far we have focused on the perturbative 



part of the axial anomaly. The perturbative QCD result ( |3.10| ) indicates that the difference 
Aqf 1 — Aqf 1 (Aq = Aq v + Aq s with the valence polaization Aq v being scheme indepen- 
dent) is induced perturbatively from hard gluons via the anomaly mechanism and its sign 
is predicted to be negative. By contrast, Aqf I (x) can be regarded as an intrinsic sea-quark 
spin density produced nonperturbatively. The well-known solution to the Z7a(1) problem 
in QCD involves two important ingredients: the QCD anomaly and the QCD vacuum with 
a nontrivial topological structure, namely the ^-vacuum constructed from instantons which 
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are nonperturbative gluon configurations. Since the instanton-induced interactions can flip 
quark helicity, in analog to the baryon-number nonconservation induced by the 't Hooft mech- 
anism, the quark-antiquark pair created from the QCD vacuum via instantons can have a 
net helicity. It has been suggested that this mechanism of quark helicity nonconservation 
provides a natural and nonperturbative way of generating negative sea-quark polarization 

In retrospect, the dispute among the anomalous gluon and sea-quark explanations of the 
suppression of I\ mostly before 1996 is considerably unfortunate and annoying since the fact 
that g±(x) is independent of the definition of the quark spin density and hence the choice of 
the factorization scheme due to the axial-anomaly ambiguity is presumably well known to 
all the practitioners in the field, especially to those QCD experts working in the area. 

3.2 Lattice calculation of the proton spin content 

The present lattice calculation is starting to shed light on the proton spin contents. An evalu- 
ation of Aqci, the gauge-invariant quark spin component defined by Agci = (p, s|g7 M 75g|p, s)s^, 
involves a disconnected insertion in addition to the connected insertion (see Fig. 2; the in- 
finitely many possible gluon lines and additional quark loops are implicit). The sea-quark 
spin contribution comes from the disconnected insertion. 



There are four lattice calculations done by Dong et al. H8], Fukugita et al. Gockeler 



et al. |50, |ST1| in the quenched approximation and Giisken et al. E2J in full lattice QCD. Note 



that the disconnected contribution is not evaluated in ^] . From Table III it is clear that 
the lattice results for g° A and As are in agreement with experiment, while the full lattice QCD 
calculations for g A , Au and Ad are too small compared to experiment. In particular, there is 
a 30% discrepancy between the lattice QCD estimate of g\ |5^| and the experimental value. 
This points to the presence of sizeable higher order or even nonperturbative contributions 
to the renormalized factor Z A S on the non-singlet current. 




connected disconnected 
Figure 2: Connected and disconnected insertions arising from the flavor-singlet axial- vector current. 

As for the chiral-invariant quantity AEqi, it involves the matrix element of J§ in light- 
front gauge where J M 5 is an anomaly-free singlet axial vector current and hence sizeable gauge 
configurations are needed in lattice calculations for AEci- Nevertheless, it is conceivable to 
have lattice results for AG and Agci soon in the near future. The lattice calculation of the 
quark total angular momentum was also available recently, see Sec. 3.4. 
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Table III. Quark spin contents of the proton from lattice calculations. Experimental results are 
taken from (2.1) for the first moments of quark spin distributions and Table II for the second 
moments^ A^q v (= Jq xAq v {x)dx) of valence quark spin densities. 





Dong et al. 


Fukugita et al. 


Gockeler et al.* 


Giisken et al. 


Expt 


9°A 


0.25(12) 


0.18(10) 




0.20(12) 


0.25(4) 


9a 


1.20(10) 


0.985(25) 


1.24(10) 


0.907(20) 


1.2670(35) 


9\ 


0.61(13) 






0.484(18) 


0.585(25) 


Au 


0.79(11) 


0.638(54) 


0.84(5) 


0.62(7) 


0.81(1) 


Ad 


-0.42(11) 


-0.347(46) 


-0.24(2) 


-0.29(6) 


-0.45(1) 


As 


-0.12(1) 


-0.109(30) 




-0.12(7) 


-0.11(1) 


A^u v 






0.198(8) 




0.169(22) |H 
0.13(1) H 


A^d v 






-0.0477(33) 




-0.055(29) H] 
-0.02(3) [f2|] 



* Since the disconnected contributions are not calculated in |5(j, |5l[], Aq receives contributions only from the 
valence component. The lattice results for A^q v , which are presumably scale independent, are obtained in 
[H at the scale fi 2 = 4GeV 2 . 



3.3 NLO analysis of polarization data 



The experimental data of gi(x, Q 2 ) taken at different x-bin correspond to different ranges of 
Q 2 ; that is, Q 2 of the data is x-bin dependent. Hence, it is desirable to evolve the data to a 
common value of Q 2 in order to determine the moments of g\ and test the Bjorken sum rule. 
Because of the availability of the two- loop polarized splitting functions APff (x) |53|, [ 



it 



became possible to embark on a full next-to-leading order (NLO) analysis of the experimental 
data of polarized structure functions by taking into account the measured x dependence of 
Q 2 at each x bin. The NLO analyses have been performed in the MS scheme (a family of the 
GI scheme) [||, ||, |7|> [17], |j| g|, ||, [yj 0, H, ||, |, ||, the Adler-Bardeen (AB) scheme 

Of course, physical quantities such as the 



yj, |65l 0, H, |6|, m and the CI scheme |7|, [ 
polarized structure function g\(x) are independent of choice of the factorization convention. 
Physically, the spin-dependent valence quark and gluon distributions should be the same in 
all factorization schemes. 

The Q 2 dependence of parton spin densities is determined by the spin-dependent DGLAP 
equations: 



d 



Aq m (x,t) 



dt 

d_ ( Aq s (x,t) 
dt \AG(x,t) 



^-AP™(x)®Aq m (x,t), 

a s {t) (APS(X) 2n f AP qG (x) 
2tt \AP Gq (x) AP GG {x) 



( Ag s (x,t) 
\AG(x,t) 



(3.20) 



with t = HQ 2 /A 2 



QCD ' 



AqNs(x) = Agi(x) - Aqj(x), Aq s (x) = ^Ag^x), 

i 



(3.21) 
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and 



AP^x) = APf(x) + ^AP«(x) + • • • . (3.22) 
The spin-dependent anomalous dimensions are defined as 

A 7 ™ = f Q \l>M)r" <l-r = A 7 if n + ^AT^'" + " " " • (3-23) 
To the NLO, AP« and AP$ were calculated in the MS scheme by Zijlstra and van Neerven 



||53|| . However, the other two polarized splitting functions APq^ and AP^g were not available 
until 1995 fl54|. The recent analysis M shows that the NLO \ 2 is significantly smaller than 



that of LO, indicating the necessity of NLO fit of data in practice. 

Since the unpolarized parton densities are mostly parameterized in the MS scheme and the 
two-loop splitting functions are available in the same scheme, it is quite natural to perform 
the NLO analysis in the MS scheme. In principle one can also work in any other factorization 
scheme. The splitting functions AP^ in the CI scheme, for example, is obtained by applying 
Eq. (|3~T0D to the spin-dependent DGLAP equation, see fl44| , f42f . It is worth accentuating 



again that though it is perfectly all right to analyze the data in the AB scheme, one has 
to keep in mind that not all hard effects are absorbed in the gluonic coefficient functions in 
such a scheme. 

The sea-quark and gluon spin distributions cannot be separately determined from current 
experimental data. In other words, while the shapes of the spin-dependent valence quark 
distributions are fairly constrained by the data, sea-quark and gluon spin densities are almost 
completely undetermined. In particular, AG is rather weakly constrained by the data. This 
is understandable because the gluon polarization contributes to g\ only at the NLO level 
through the convolution AGg <8> AG. In principle, measurements of scaling violation in 
gi(x, Q 2 ) via, for example, the derivative of gi(x, Q 2 ) with respect to Q 2 , in next-generation 
experiments will allow an estimate of the gluon spin density and the overall size of gluon 
polarization. Of course, the data should be sufficiently accurate in order to study the gluon 
spin density. Meanwhile, it is even more important to probe AG (a;) independently in those 
hadron-hadron collision processes where gluons play a dominant role. 

Though the polarized structure function is factorization scheme independent, it is impor- 
tant to perform NLO analyses in different schemes to test the reliability and consistency of 



the theory. It is found in |L9], |58|, |y] that the polarized parton distributions obtained in MS, 
AB and CI schemes agree well for the non-singlet spin densities, while the first moment of 
AG{x) obtained in the AB or CI scheme is different from that in the MS scheme, reflecting 
that the present data can hardly constrain the gluon spin distribution. Typically, the ex- 
tracted value of the gluon spin at Q 2 = 1 GeV 2 lies in the range < AG < 2 . However, the 



recent analysis |64| shows that the LO fit cannot decide on the sign of AG, while the NLO 
analysis yields a negative first moment of the gluon density^ This illustrates again that it is 
difficult to pin down the gluon spin distribution from present polarized DIS data. Note that 

4 It has been found by Jaffe |36| that the gluon spin component is negative AG ~ —0.4 in the MIT bag 
model and even more negative in the non-relativistic quark model. However, it was explained in f6^| that 
the negative AG obtained by Jaffe is a consequence of neglecting some self-interaction effects. 
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a recent NLO analysis in shows that the polarized strange quark density is significantly 
different from zero independently of the factorization schemes used in the analysis: 



As m = -0.102 ± 0.012, Asci = -0.064 ±0.010, As AB = -0.058 ± 0.012 , (3.24) 

at Q 2 = 1 GeV 2 . Note that the sea polarization As is scheme dependent: 

ot 

As C i,ab = As GI + -^AG. (3.25) 

The presence of the sea polarization in the CI or AB scheme implies that the gluon polar- 
ization is not at its maximal value given by AG ~ — (27t/qi s )Asgi- 

3.4 Orbital angular momentum 

The orbital angular momentum plays a role in the proton spin structure. For example, the 
growth of large AG with Q 2 is balanced by the large negative orbital angular momentum of 
the quark-gluon pair. It is also known that the reduction of the total spin component AE 
due to the presence of the quark transverse momentum in the lower component of the Dirac 
spinor is traded with the quark orbital angular momentum. 
In the proton spin sum rule: 

\ = J q + Jg, (3-26) 

the total angular momenta J q and Jq of quarks and gluons respectively are gauge invariant. 
However, the decomposition into spin and orbital components, J q = | AS + L q and Jg = 
AG+Lq, is gauge dependent. This leads to difficulties in defining the partonic spin densities: 
There exist no local gauge invariant operators that could represent the densities of gluon 
spin and the orbital angular momenta of quarks and gluons. It is known that the spin and 



orbital angular momenta of quarks and gluons appearing in the decomposition [[34 



J z — J n + Jn — S* + L z n + S%< + L) 



d 3 x 



G ^G 

1 



-V^7 3 7 5 ^ + ip\x x id) 3 ip + {Ex A) 3 - E k {x x d) 3 A k 
2 



(3.27) 



obtained in the the light-cone gauge A + = and in the infinite momentum frame, are sep- 
arately gauge variant except for the quark spin operator S q . However, the gluon spin and 
its distribution, for example, are physical, gauge invariant quantities and can be measured 
experimentally. The point is that AG can be expressed as the matrix element of a string-like 
nonlocal gauge- invariant operator Oq p4| , |68| . As stressed in |59|, what one measures exper- 



imentally is the matrix element of the gauge-invariant operator Og- But in the light-cone 
gauge, the above operator reduces to the gluon spin operator S G = J d 3 x(E x A) 3 ; that is, the 
gauge-invariant extension of the gluon spin operator in the light-cone gauge is measurable. 
Likewise, a gauge-invariant operator that reduces to the quark orbital angular momentum 
in the light-cone gauge has been discussed in [7(J (see however a different discussion in |7T| ). 



In principle, the total angular momenta of quarks and gluons, J q and Jq respectively, 
can be measured in deeply virtual Compton scattering in a special kinematic region where 
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single quark scattering dominates [ff^] . It has also been suggested that the orbital angular 
momentum might be deduced from an azimuthal asymmetry in hadron production with a 
transversely polarized target. 

The evolution of the quark and gluon orbital angular momenta was first discussed by 



equation in the light-cone gauge: 



Ratcliffe |[73[ . Ji, Tang and Hoodbhoy |74| have derived a complete leading- log evolution 

(3.28) 






tJ \ ag J 



with the solutions 



LJQ 2 ) = --AE + - 3Uf + f(Q 2 ) [ LJQl) 
gvv ) 2 2 16 + 3n/ ) y <A^o, 

L G {Q 2 ) = -AG{Q 2 ) + I- + f(Q 2 ) (l g {QI) + AG(Q 2 ) 



2 16 + 3n 



/ 




where 



\nQ 2 /A 2 

m > = i^C 1 (3 30) 

and AS is Q 2 independent to the leading-log approximation. We see that the growth of 
AG with Q 2 is compensated by the gluon orbital angular momentum, which also increases 
like \\iQ 2 but with opposite sign. The solution ( |3.29D has an interesting implication in the 
asymptotic limit Q 2 — > oo, namely 



J q (Q 2 ) = -AS + L q (Q 



2 ' 2 16 + 3n/ ' 

1 1 fi 

J G (Q 2 )=AG(Q 2 ) + L G (Q 2 )^ 2 16T3^- (3 ' 31) 

Thus, history repeats herself: The partition of the nucleon spin between quarks and gluons 
follows the well-known partition of the nucleon momentum. Taking rif = 6, we see that 
Jq '■ Jg — 0.53 : 0.47 . If the evolution of J q and Jq is very slow, which is empirically known 
to be true for the momentum sum rule that half of the proton's momentum is carried by 
gluons even at a moderate Q 2 , then AS ~ 0.25 at Q 2 = 10 GeV 2 implies that L q ~ 0.13 at 
the same Q 2 , recalling that the quark orbital angular momentum is expected to be of order 
0.20 in the relativistic quark model. 

Recently the quark orbital angular momentum of the nucleon was calculated from lattice 



QCD by considering the form factor of the quark energy-momentum tensor T^ v |75| . The 
total angular momentum of the quarks is found to be 

J q = 0.30 ±0.07. (3.32) 

That is about 60% of the proton spin is attributable to the quarks. Hence, the quark orbital 
angular momentum is L q = 0.17 ±0.06 when combining with the previous quark spin content 
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^AE = 0.13 ± 0.06 Therefore, about 25% of the proton spin originates from the quark 



spin and about 35% comes from the quark orbital angular momentum. 

It must be stressed that J q should be factorization scheme independent. This means 
that a replacement of AE G i by AEci in the spin sum rule ( |3.26| ) requires that the difference 



ASci — ASqi = (nfa s /2ir)AG be compensated by a counterpart in the gluon orbital angular 
momentum: 



rci = r Gi _ !^£ AG . (3.33) 



It is interesting to note that if AG is of order 2.5 , one will have AEci ~ 0.58 and L^ 1 ~ . 
In other words, while AEci is close to the relativistic quark model value of AS, L^ 1 deviates 
farther from the quark model expectation L' q ~ 0.20 (see Sec. I). 

4 Conclusions 

The spin sum rule of the proton in the infinite momentum frame reads 

I = J q + Jg = ^AE + L q + AG + Lq. (4.1) 

The quark spin can be inferred from polarized DIS measurements of gi(x) and its first 
moment. Due to the ambiguity arising from the axial anomaly, the definition of the sea 
polarization Aq s and hence Aq = Aq v + Aq s is k± factorization dependent, but J q , gi(x) 
and Ti are not. The only spin content which is for sure at present is the observed value 
AS ~ 0.20 — 0.30 in the GI scheme (e.g. the MS scheme). The recent lattice calculation 
yields J q = 0.30 ± 0.07. Therefore, we have Lf ~ 0.10 ± 0.06 for AS GI ~ 0.25, to be 
compared with the quark model prediction L' ~ 0.20 . The values of AS and L q in the CI 
scheme (e.g. the improved QCD parton model) or the AB scheme depend on the gluon spin. 
Since L^ 1 — L^ 1 = — nja s AG/(47r), it is clear that if the gluon polarization is positive and 
large enough, then L q 1 will deviate even farther from the quark picture although AEci can 
be made to be close to the constituent relativistic quark model. In the asymptotic limit, 
J q (oo) = |AE(oo) + L q (oo) ~ j and Jg(oo) = AG(oo) + Lg(oo) ~ \. The recent lattice 
result J q = 0.30 ±0.07 at a moderate Q 2 seems to suggest that the evolution of J q and hence 
Jg is slow enough. 
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